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ABSTRACT 

We investigate the changes in polarization position angle in radiation from 
pulsar A around the eclipse in the Double Pulsar system PSR J0737-3039A/B at 
the 20 cm and 50 cm wavelengths using the Parkes 64-m telescope. The changes 
are ~ 2a during and shortly after the eclipse at 20 cm but less significant at 50 cm. 
We show that the changes in position angle during the eclipse can be modelled by 
differential synchrotron absorption in the eclipse regions. Position angle changes 
after the eclipse are interpreted as Faraday rotation in the magnetotail of pul- 
sar B. Implied charge densities are consistent with the Goldreich- Julian density, 
suggesting that the particle energies in the magnetotail are mildly relativistic. 



Subject headings: binaries: eclipsing — polarization — pulsars: individual (PSR 
J0737-3039A, PSR J0737-3039B) 
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Introduction 



PSR J0737-3039A/B, also known as the Double Pulsar system, is a compact binary 
sys tem consisting of tw o ne utron stars, both r adiating; at radio frequencies^. Discovered 
by iBurgay et al.l (120031 ) and iLyne et al.l (120041 ) in a survey at 20 cm wavelength using a 
multibeam receiver on the 64-m radio telescope at Parkes in New South Wales, Australia, the 
two pulsars, referred to as pulsar A and pulsar B respectively, orbit each other with a period 



of about 2.45 hours in an almost circular trajectory with eccentricity 0.088 (IBurgay et al. 



20031 ; Kramer et al.ll2006l ). The inclination angle between the line of sight and orbital plane 
normal is estimated to be 87.7°. The spin periods for pulsar A and B are Pa = 22.7 ms 
and Pb = 2.77 s, and their surface dipole magnetic field strengths are E>a = 6.3 x 10 9 G and 
B B = 1.2 x 10 12 G. The spin-down powers for pulsar A and B are Ea ~ 5.9 x 10 33 ■ 



erg s 



and Eb 



1.6 x 10 30 ergs 1 , 



giving Ea/Eb ~ 3500. 



Because of the strong relativistic wind from pulsar A, the magnetosphere of pulsar B is 
truncated compared to that of an isolated pulsar. The relativistic particle flux in the wind 
penetrates into the light cylinder of B resulting in the formation of a bow shock that strongly 
perturbs the dipolar field structure by compressing B's magnetosphere on the "daytime" side 
(the side facing A), and stretching it on the "nightime" side into a magnet otail. From studie s 
on simulating the interaction between A's wind and B's magnetosphere (lArons et al.l 120041 ). 
the magnetotail probably reaches a distance of several semi-major axes of pulsar B's orbit 
(the semi-major axis is 4.5 x 10 8 m). Figured] illustrates the orbital arrangements for the two 
pulsars and the magnetotail. The wind from pulsar A blows at pulsar B's magnetosphere 
causing B's magnetospheric material to flow in the radial direction, but because of the orbital 
motion, the magnetotail is curved behind pulsar B. 

The fact that the orbit is almost edge-on to the line of sight makes an eclipse observable 
each orbit when pulsar A moves behind pulsar B, as shown in Figure [TJ The eclipse occurs 
at orbital longitud e of ~ 90° measured from the ascending node and lasts for ~ 30 seconds 
(ILyne et al.l 120041 ) . It was shown that, during an eclipse, the pulsed fl ux density of A is 



modu lated at the rotational period of pulsar B and its second harmonic ([McLaughlin et al. 



20041 ) . Synchrotron absorption by the highly relativistic plas ma in the closed field-line region s 
of pulsar B has been proposed as the eclipse mecha nism (ILyutikov and Thompson! 120051 ). 
Modelling of these variations by iBreton et al.l (120081 ) gave best-fit values of ~ 130° for the 
angle between rotation axis and the orbital normal, and ~ 70° for the inclination angle 
between the rotation and magnetic axes. Hence pulsar A's radiation passes through pulsar 



1 PSR J0737-3039 B is currently unde tectable at radio frequencies, probably because of relativistic preces- 
sion of its spin axis ( Perera et al. 2010l ). 
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Fig. 1. — Orbital arrangements of the two pulsars and the magnetotail behind pulsar B (not 
to scale). The black and grey curves represent the orbits for pulsar A and B respectively and 
both pulsars orbit in the clockwise directio n. The shade d areas around pulsar B represent 
the closed field-line regions (adapted from iBreton et al.l (120081 ) and online material), and 
bj and /i are the angular velocity and magnetic moment vectors respectively. The orbital 
normal is directed out of the page. As the two pulsars orbit, the radiation we observe from 
pulsar A first encounters the leading edge of the shock front, enters the magnetosphere of 
pulsar B, leaves the magnetosphere on the other side and finally crosses the magnetotail. 
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B's truncated magnetosphere and interacts with the plasma during an eclipse, then, as the 
two pulsars advance in orbit away from the superior conjunction, pulsar A's radiation begins 
to traverse the magnetotail. 

The radiation from pulsar A that penetrates the magnetosphere of pulsar B during and 
shortly after an eclipse serves as a probe into the plasma content of B's magnetosphere and 
magnetotail. Study of the emergent radiation may reveal the largely unknown properties 
of pulsar magnetospheres and possible physical processes within them. In this paper, we 
investigate the changes in polarization position angle in radiation from pulsar A at 20 cm 
and 50 cm wavelengths during and around the eclipses. 



2. Observation and data analysis 



Our observations of the Double Pulsar use the Parkes 64-m radio telescope. The data 
of interest were taken in two frequency bands centred at 1369 MHz (20 cm) from September 
2007 to October 2011, and at 730 MHz (50 cm) from February 2010 to November 2011. 
The 20 cm observations wer e made using the center beam of the 20 cm Multibeam receiver 
( jStaveley-Smith et al.l Il996f) and the 50 cm observations used the 10 cm/50 cm dual-band 



receiver f lGranet et al.l 120051 ). The Parkes digital filterbank systems, PDFB3 and PDFB4, 
were used to split the 20 cm and 50 cm data into 1024 and 512 channels respectively, form 
the polarization products and fold the data at the topocentric pulse period. All observations 
had a sub-integration time of 30 s. We checked the presence of an eclipse in a data file 
by checking the coverage of orbital phase of 90°. Once we had identified all relevant data 
files, regions lOmin before an d after the eclipse were retained for further analysis using 
the psrchive software suite ( iHotan et al.l 120041 ). Spectral band edges were deleted and 
interference was removed from all files. Files that had significant interference at or close to 
the eclipse region were eliminated from our list. Short observations of a linearly polarized 
broad-band calibration signal preceding each pulsar observation were used to calibrate the 
instrumental phase and gain. For the 20 cm system, the effects of cross-coupling between the 
orthogonal feed probes were removed as part of the calibration pro cess. Data were summe d 



in frequency using the corrected rotation measure, + 112.3radm" 2 (Manchester et al.ll2010h . 
We obtained 38 files at 20 cm and 13 at 50 cm with good quality data. 

To improve the signal-to-noise ratio, we summed the files in each band, binning the 
results in orbital phase. The phase bins were ~ 0.7° or ~ 18 s in width. Sub-integration 
data were assigned to the two adjacent orbital phase bins with weights proportional to 1 — x 
and x for the leading and trailing bins respectively, where x is the interval between the 
leading phase bin centre and the sub-integration center time, normalized by the phase bin 
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width. Finally, Stokes parameters were summed across the main pulse and interpulse to form 
the average Stokes parameters for each orbital phase bin. This summing of data across the 
two pulses is possible because the strongly linearly- polarized parts of both pulse components 
have essentially constant and equal position angle (IDemorest et al.l 12004 ) . 



3. The results 

Measurements of the mean position angles and Stokes I (total intensity) for 20 cm and 
50 cm wavelengths as a function of orbital phase are presented in Figure [2J Uncertain- 
ties in the position angle f or each orbital phase were computed using a = 28.65° (aj/L c ) 



( jLorimer and Kramerll2005l ). where 07 is the standard deviation of the total intensity, and 
L c is the linear polarization corrected for noise bias. The Stokes I curves, which are scaled 
to fit into the plot, clearly show a large drop in intensity representing the eclipse region. 

We consistently detected changes in position angle during and shortly after an eclipse 
at both wavelengths. Inspection of the 20 cm curve in the eclipse region shows that the 
position angle decreases during the ingress phase with a change of approximately 6° ± 3° 
relative to the mean baseline level. It then rises to its maximum value of approximately 
3° ± 2° above the baseline at the end of the eclipse. It remains high for the following phase 
bin, then drops back to the baseline value. There are also changes in polarization angle 
in the corresponding regions during the eclipse at 50 cm, but the uncertainties are larger. 
Nevertheless, the position angle changes shortly after the eclipse is consistent with those at 
20 cm. 

We confirmed our results by randomly partitioning the 20 cm files into two subsets 
as shown in Figure The changes in position angle during the eclipse are essentially 
consistent between the two subsets, i.e., a decrease during ingress followed by an increase 
during egress. Partitioning the dataset into either two, three or four subsets, either randomly 
or chronologically, produced similar results. 



4. Mechanisms to change the polarization position angle 

4.1. Synchrotron absorption 

Synchrotron absorption, like synchrotron emission, may be separated into two linearly 
polarized components, with the more strongly emitted and absorbed component having its 
electric vector perpendicular to the magnetic field. The polarization-averaged synchrotron 
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Fig. 2. — Plot of polarization position angle in degrees (upper curves) and Stokes I (total 
intensity) parameter (lower curves) plotted against orbital phase in degrees for 50 cm (solid 
curve) and 20 cm (dashed curve) at ~ 18 s time resolution. The curve for position angle at 
50 cm is offset by 3° for clarity. The Stokes I curves clearly show the eclipse centred near 
90°. The horizontal lines represent average baseline values. 
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Fig. 3. — Position angles at 20 cm for two subsets of the total data set plotted against orbital 
phase. 



-8 - 



absorption coefficient for a rela tivistic thermal distribution of particles, n(e) oc e £ / fcsT , is 
(jLyutikov and Thompson! 120051 ) . 



a,, 



Air' 2 en s 
3 7 / 3 B sin k 



-i 5/3 



v-q e sm k 



u 



for u <C Up ~ 4(kBT e /m e c 2 ) 2 UB, e , where n s , B, m e , c, &b, T, and u are the density of 
the absorbing particles, magnetic field strength, electron mass, speed of light, Boltzmann 
constant, plasma temperature and radiation frequency respectively, k is the angle between 
the magnetic field and propagation vector of the radiation and UB, e = eB /2irm e c is the 
cyclotron frequency. The absorption coefficients for polarization perpendicular and parallel 
to the magnetic field are aj; = 4a u /3 and al = 2a u /3, and so a^/al = 2. The optical 
depth is t v = a u L, where L is the effective path length through the absorbing medium. 
Hence, synchrotron absorption preferentially removes the linear polarization perpendicular 
to the magnetic field. As a result, the position angle of the transmitted radiation changes 
with increasing optical depth, with the position angle approaching the angle of the projected 
magnetic field. If Xo an d Xu are the incident and emergent polarization angles relative to 
the projected magnetic field direction, then, 



1^ = I„smxv = Jo^expf-r^) = I 0)i , sin^o exp(-4r I/ /3) 
/J] = i^cosx^ = /(S^exp(-rJ) = I 0) „ cos% exp(-2r v /3). 



(2) 



where J -1 '" represents intensity perpendicular and parallel to the projected magnetic field 
respectively, and quantities with subscript zero represent incident values at a particular 
frequency, u. The change in position angle is given by A^ u = (xu — Xo)- 



4.2. Faraday rotation 

The decomposition of linearly polarized radiation into two circular modes that propagate 
at slightly different phase velocities results in Faraday rotation of the plane of the received 
linear polarization. Integrating over the full path gives the change in position angle, 

Atp ps 0.81 A 2 / n e B cos 9 ds = 0.81L (n e B) A 2 = RM A 2 , (3) 
Jo 

where n e , B and L are the electron density in cm -3 , the average magnetic field in microgauss, 
and the pathlength in parsecs respectively, and cos 9 is the angle between the magnetic field 
and the ray path, RM is the rotation measure. The notation (...) implies an average along 
the ray path. 
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Pulsar magnetospheres are thought to be populated with highly relativistic plasma. 
If there is no thermal plasma present, as some models imply, the modes would be nearly 
linearly polarized, and conventional Faraday rotation would not occur. The natural wave 
modes would also be linear if the density of the electrons and positrons were the same. The 
circular component needed can arise from the net charge, associated with the Goldreich- 
Julian density, provided the particle energies are not highly relativistic. Any measured 
Faraday rotation would have the interesting implication that the dominant contribution to 
the wave dispersion is from electrons that are no more than mildly relativistic. 



5. Discussion 

5.1. Changes in position angle during the eclipse 

Our measurements of the position angle change of A's pulse emission integrate over the 
rotation of B's magnetosphere and the resulting eclipse modulation, with our 18 s averaging 
time covering approximately 6.5 revolutions of B. The direction of the differential synchrotron 
absorption relative to the incoming radiation will vary during B's rotation and along the 
path through the magnetosphere. It will also change as the line of sight to A moves through 
B's magnetosphere because of the orbital motion. Our present results do not justify a 
full computation of this complex variation. However, we can approximate the effect by 
recognising that the closed field lines perpendicular to the magnetic axis are symmetrical 
but of opposite sign in the two hemispheres. As the star rotates, they will therefore average 
to zero leaving a net magnetic field along the magnetic axis. Furthermore, since the magnetic 
and rotation axes are not aligned, the field perpendicular to the rotation axis will also cancel 
leaving a net magnetic field along the rotation axis of B. During the four years covered by 
our observations, B's spin axis will precess by about 20° but this will have little effect of the 
observed absorption which is dominated by averaging over B's rotation 



The peaks and dips in flux density from pulsar A in the eclipse light curve ([McLaughlin et al 



20041 ) imply that the optical depth is, respectively, either ~ when the line of sight misses 
B's magnetopsheric torus, or > 1, which corresponds to the absorbing closed- field region 
that obstructs our view of pulsar A resulting in a large decline in flux density. The differ- 
ential changes in /J - and jj would be insignificant during most of the eclipse region where 
r > 1 because the intensities are low, except near the boundaries of the eclipse where 
the optical depth r ~ 1, which characterizes the applicability of Equation (j5J). For the 
synchrotron absorption model (Equation [I]), the density of the absorbing plasma can be 
esti mated by n 9 = Mng . j, whe re M = 10 6 is the pair multiplicity parameter, and uqj is 



the iGoldreich and Julian! ( 119691 ) charge density. The energy balance between A's relativistic 
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wind and B's magnetic pressure gives a field strength of 7 G at t he eclipse bounary, which 
is est imated to be near r\ rA \ « r^/10 ~ 1.5 x 10 9 cm from pulsar B (ILyutikov and Thompson 



20051 ). and so uqj ~ 0.2 cm -3 , giving n s ~ 2 x 10 5 cm -3 . The cyclotron frequency, vs,e is 
2 x 10 7 Hz, and the relativistic electrons are assumed to be thermally distributed at temper- 
ature [ksT /m e c 2 ) = 10. Assuming B's magnetic field is randomly oriented and averaging 
between 0° and 90° gives ((sulk) 2 / 3 ) »s 2/tt, and with L = 7"b a i, we obtain (r 2 o) ~ 1 and 
(r 50 ) ~ 3, where the number in subscript represents the wavelength. Using these values and 
our data for the total intensity, and assuming xo — 45°, Equation (J2]) gives an order of mag- 
nitude change in polarization or position angle of 10°, consistent with our measurements. We 
therefore, conclude that differential synchrotron absorption in the closed field-line regions 
can plausibly account for the observed position angle changes during the eclipse. 

Combining Equation ([2]) and rearranging gives tan Xu = tan (xo+Aip u ) = tan xo exp(—2r u /3). 
In principle, with higher time-resolution observations of the eclipse, one could use this equa- 
tion to determine the angle of the rotation axis of pulsar B projected on the sky, ips,n- At 
times when r v is neither very large nor very small, both r„ and Aip u may be measurable. 
Since Xo — ipo ~ ipB,n and Aip u = Xu — Xo = tyv ~ ^o> the above equation can be rewritten 
as tan(^„ — ipB,n) = tan(^ ~ ipB,n) exp(— 2t u /3). With multiple measurements for each of 
tp u and t v , and since ipo is known, ipB,n may be determined. 



5.2. Changes in position angle due to Faraday rotation 



After the eclipse, the synchrotron absorption effect will diminish and Faraday rotation 
should start to manifest. During this time, pulsar A's radiation goes through regions located 
immediately beyond the closed field-line boundaries in the wind zone where the oppositely 
directed magnetic fields in the two lobes from the magnetotail each connect to one of the 
polar caps. Unlike the g eomagnetotail where field structure is symmetric in these regions 
(IPilipp and Morfilllll978l ). the large inclination angle results in an asymmetric field structure 
between the two lobes in these regions and an alternating field stru cture as the pulsar rotates 
within the confining magnetosheath (jSpitkovsky and Aronsll2004l ). With the misalignment 
between the rotation axis and the orbital normal, a net magnetic field in the line of sight will 
result in Faraday rotation in radiation that traverses these regions of pulsar B's magnetotail. 

Faraday rotation causes changes in position angle of the same sign at both wavelengths 
shortly after an eclipse. The measured changes in position angle can be used to estimate 
(n e B)L, and a separate estimate of L and B allows n e to be estimated in the magnetotail. 
Consider mean values along the line of sight, and assuming magnetic field scales as oc r~ 2 up 
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to L/2 in the magnetotai@ gives B ~ 2 G, where L « 1.5 x 10~ 8 pc is the semi-major axis 
of pulsar B's or bit. Equation (|3l), with A^ n ~ —5°, gives n e ~ 40 cm -3 . The corotating 
charge density (jGoldreich and Julian! Il969l ) is the difference between the number densities 
of positrons and electrons, which gives ~ 20 cm" 3 at r = r^ai- Although the corotation 
requirement does not apply beyond rbai, the wind in the magnetotail rotates with the star 
implying a charge density of order the Goldreich- Julian value. This suggests that Faraday 
rotation may be responsible for the changes in position angle away from the eclipse region. As 
mentioned earlier, the complication in interpreting the results lies in the complex variations 
in the orientation of B's magnetic field structure, caused by field reversals between the two 
lobes in the magnetotail, the large inclination angle between magnetic and rotation axes of 
pulsar B and the orbital motion of the two pulsars. However, with higher time resolution, it 
is possible that changes in A's position angle could be correlated with different orientations 
of B's field structure. 



5.3. Mildly relativistic particle energies 



Faraday rotation in the magnetotail would require that the charged particles be no 
more than mildly relativistic. If the particles were highly relativistic, the natural wave 
modes would be linearly polarized and there would be no Faraday rotation. If our sugges- 
tion that significant Faraday rotation occurs is confirmed by more detailed observations, it 
would be a strong indication that the bulk of the pairs produced in a pulsar magnetopshere 
are no more than mildly relativistic. This is inconsistent with theories for t he pair creation 
that suggest much higher mean Lorentz factors of typically 7 



Ruderman and Sutherland 1975 ; Daughertv and HardineJ 



10 2 - 10 4 ( Sturrocklll971 



1982)^ S ince energy-loss mecha- 



nism s are ineffective (IRuderman and Sutherlandlll975l ; lSturnerlll995l ; iLyubarskii and Petrova 
20001). any nonrelat ivistic plasma must be created directly in the pulsar magnetosphere 



fjJessner et al.ll200ll ). Direct generation of mildly relativistic pairs is favored for magnetic 

i 4.9 x 10 13 G is the critical field strength 



fie lds with 0.02^ < B < 0.1 B cr , where B, 



flWeise and Melrose! l2002t ) 



2 Scaling of the magnetic field strength in the open field-line region and magnetotail is uncertain but not 
critical to our discussion. 
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6. Conclusions 

We have confirmed the changes in polarization position angle in radiation from pulsar 
A during and shortly after the eclipse in the Double Pulsar system at 20 cm and 50 cm 
wavelengths. The changes are ~ 2a at 20 cm and ~ la at 50 cm. The variations at 20 cm were 
confirmed by partitioning the data files into two subsets. Differential synchrotron absorption 
in the closed field-line regions can account for the observed position angle changes during 
the eclipse. Modelling of the changes in position angle shortly after the eclipse with Faraday 
rotation gives a charge density that is consistent with the Goldreich- Julian value. Since 
energy-loss mechanisms in pulsar magnetospheres are inefficient, the presence of Faraday 
rotation implies that pairs must be created directly in nonrelativistic regime. 
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